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ABSTRACT 

Observational measurements of the relationship between supermassive black holes (SMBHs) and the prop- 
erties of their host galaxies are an important method for probing theoretical hierarchical growth models. Grav- 
itational lensing is a unique mechanism for acquiring this information in systems at cosmologically significant 
redshifts. We review the calculations required to include SMBHs in two standard galactic lens models, a cored 
isothermal sphere and a broken power law. The presence of the SMBH produces two primary effects depending 
on the lens configuration, either blocking the "core" image that is usually predicted to form from a softened 
lens model, or adding an extra, highly demagnified, image to the predictions of the unaltered lens model. The 
magnitudes of these effects are very sensitive to galaxy core sizes and SMBH masses. Therefore, observations 
of these lenses would probe the properties of the inner regions of galaxies, including their SMBHs. Lensing 
cross-sections and optical depth calculations indicate that in order to fully observe these characteristic signa- 
tures, flux ratios of order 10 6 or more between the brightest and faintest images of the lens must be detectable, 
and thus the next generation of radio telescope technology offers the first opportunity for a serious observa- 
tional campaign. Core images, however, are already detectable and with additional observations their statistics 
may be used to guide future SMBH searches. 
Subject headings: black hole physics — galaxies:evolution — gravitational lensing — surveys — telescopes 



1. INTRODUCTION 

The central regions of galaxies contain a great deal of in- 
formation regarding their formation and evolution. The mass- 
density profiles of these inner regions have been recognized 
for sometime as importan t measures of cold dark matter mod- 
els (iNavarroetaU 119971 [Keeton & Madau 2001), and more 
recently, the properties of central supermassive black holes 
(SMBHs) have emerged as important probes into the history 
of hierarchical growth models. 

SMBHs are inferred compact objects with masses between 
10 6 and 10 9 solar masses. Several empirical correlations 
have been identified between SMBH masses and the prop- 
erties of their host galaxies, including correlations with the 
mass or luminosity of the galactic bulge ( Magorrian et al. 
119981 iKormendv et al.ll2000t Laor | l2001l) . with the bulge ve- 
locity dispersion (Gebhardt et al. 2000; Merritt & Ferrarese 
2001), and also wit h the light concentration of the bulge 
(Grah am etafl feOOl ). Recent analysis suggests that SMBH 
mass may additionally be correlated with the circular ve- 
locity of the ho st galaxy, even beyond the optical radius 
(Ferrarese 2002; Baes et al. 2003), and thus there may be a 
link between the SMBH mass the total mass of the galaxy. 
Because of this apparently tight link between SMBHs and 
their host galaxies, including dark matter halos, it is likely 
that the formation of both were closely related. Most 
contemporary hie rarchical growth models do in fact re- 
flect this property (Silk & Rees 1998; Kauffmann & Haehnelt 
2000; Burkert&Silk 2001; MacMillan & Henriksen 2002; 
Wvith e & Loebll2002t IVolonteri et al.ll2003t iDi Matteo et al.l 



2003 

The detection of SMBHs is difficult and traditionally has 
been achieved using stellar kinematics in the very central 
regions of galaxies. These techniques have been feasible 
only for relatively nearby galaxies (D < 100 Mpc, z < 0.02) 
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llTremaine et alJ [2002) in which the cores can be resolved. 
SMBHs are presumed to be the inactive relics of active galac- 
tic nuclei, and thus, have evolved over a long redshift his- 
tory (z ^ 6 to present). Motivated, in part, by the desire 
to observe this evolution, additional techniques for deduc- 
ing SMBH mass in more distance galaxies have been devel- 
oped. The most prominent of these is the use of gravitation- 
ally broadened emission lines in active galactic nuclei t o infer 
SMBH mass in these objects (M cLure & DunlorJ20'0ll) . 

Gravitational lensing offers a compelling mechanism to ex- 
plore the relationship between the properties of SMBHs and 
their host galaxies at intermediate redshifts, extending to z — 1 
or greater. Examination of strong galaxy-galaxy (or galaxy- 
quasar) gravitational lenses under proper circumstances may 
yield the mass of the SMBH, information on the inner mass- 
density profile of the host galaxy, and the total mass of the 
galaxy simultaneously, greatly simplifying existing observa- 
tional hurdles. 

Relatively minimal consideration has been given to the 
gravit ational lensing properties of SMBHs in galaxies, how- 
ever. Alexander & Loeb (2001) have considered the SMBH 
in the Milky Way, discussing the influence of stellar neigh- 
bors on its lensing properties, as well as the effects of its 
gravitational shear on the microlensing of individual stars. 
For galaxy-galaxy lensing, only Mao et al. (2001) and Chen 
(2003) have made general analyses. The former demon- 
strate calculations for adding a SMBH to a standard galac- 
tic lensing potential and the latter s hows how SM BHs can 
affect image separations in lenses. Keetorj ll2003l) consid- 
ers a sample population of 73 mock lenses and calculates 
lensing statistics including SMBHs effects. He is inter- 
ested in observing core images, which contain information 
about the inner mass-density profiles of galaxies, and finds 
that SMBHs should reduce the probability of observing the 
faintest of these images, but otherwise not generally affect 
core images. Using observational data, Rusin & Ma (2001) 
show that explicitly including SMBHs in models weakens 
constraints on mass profiles in the inner regions of galax- 



ies and Winn et al. (2004) have achieved an upper bound 
of 2 x 10 8 solar mas ses for any SMBH in the l ens s ystem, 
PMN J1632-0033. Wambsganss & PaczvnskJ dl992l) dis- 
cussed using VLBI observations to detect SMBHs in galaxy 
halos, and Wilkinson et al. (2001) and Augusto & Wilkinson 
(2001) placed limits on the abundance of supermassive com- 
pact objects in the universe. 

In this paper we expand upon these initial efforts and pro- 
vide a general context for interpreting SMBH effects in lens- 
ing studies. In particular, we will elaborate on the discus- 
sions of Keeton (2003) and Mao et al. (2001) regarding the 
impact of SMBHs on observing core images and also per- 
form calculations to determine the likelihood of directly ob- 
serving SMBH lensing signatures. We begin in Section|2]by 
presenting two common galactic lens models and perform- 
ing the calculations to describe the addition of a SMBH to 
each. From these modified models we establish criteria which 
a lens-source system must meet in order for SMBH effects 
to be detectable. In Section [3] we use these criteria to com- 
pute cross-sections for the production of observable lensing 
events. Using the observed local relation between SMBH 
mass and bulge velocity dispersion, these cross-sections are 
convolved with the velocity dispersion distribution of galaxies 
to establish estimated optical depths to lensing. We conclude, 
in Sections |4] and |5j by summarizing the lensing signatures 
of SMBHs and discussing their effects on other lensing stud- 
ies, such as the search for core images, and also how SMBH 
signatures themselves might be observed with existing and 
proposed future facilities. 

2. ANALYTICAL TREATMENT OF LENS MODELS 

2.1. Lensing by a Single Galaxy with a SMBH 

In order to determine the statistics of lensing by a super- 
massive black hole in a galaxy core, we must first address the 
mechanics of the lensing. In this section we review the prop- 
erties of combined galaxy and black hole lensing, emphasiz- 
ing those aspects that are important for the determination of 
cross-sections for lensing statistics. The SMBH is modeled 
as a point mass and we consider two models for the poten- 
tial of a lensing galaxy, a cored isothermal sphere and a bro- 
ken power law. It has been recognized for some time that 
the observable surface brightnesses in the central regions of 
galaxies do not actually approach zero slope as implied by 
the canonical cored isothermal potential. This has been the 
motivation for various broken power law parameterizations 
for surface brightness, such as the "Nuker Law" (Lauer et al. 
1995; Bvun et alJl996tlFaber et alll997UZhaoll997l) . While 
surface brightness has indeed been shown to be well repre- 
sented by a broken power law, the actual mass-density distri- 
bution remains largely unknown and hence we will treat both 
models in parallel for this paper. 

2.1.1. Point Mass and Cored Isothermal Sphere (PM+ CIS) 

We begin by considering a point mass potential added to 
a cored isothermal sphere potential. A similar treatment has 
already been performed by Mao et al. (2001) starting from a 
cored isothermal sphere mass-density. The distinction is sub- 
tle and the resulting behavior is quite similar. The advantage 
to the approach followed here is that the characteristic equa- 
tions are straightforward. 

Assuming the lensing occurs in the weak fiel d limit, the 
standard one-dimensional lens equation is ( Scheider et al. 
1992, Section 5.1) 

/? = 0-Vtf, (1) 



where 9 is the position of the image, j3 is the position of the 
source as measured from the center of the lens, \& is the two- 
dimensional projected potential of the lens, and V^ repre- 
sents the deflection angle. For an isolated point mass, 



y = 6 l PM lnt 



(2) 



and for an isolated galaxy with a cored isothermal lensing po- 
tential, 



* = e sls Je 2 +9 2 c 



(3) 



(Naravan & Bartelmann 1996, Table 2). The core radius of 
the galaxy potential is denoted by 9c and the Einstein radii 
of the black hole and galaxy are denoted by 6*pm and 6*sis, 
respectively: 
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The velocity dispersion of the galaxy is er, the mass of the 
black hole is Mbh, and the distances from observer to lens, 
lens to source, and observer to source are £>ol, ^ls, and Dos, 
respectively. 

For the weak field limit, lensing potentials add linearly, thus 
a reasonable projected potential for lensing by a galaxy with 
an SMBH is obtained by simply combining these two poten- 
tials to give 



^ = 9^\ 9 2 + 9 2 + 9 2 lntf, 
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For circularly symmetric lens potentials the magnification at 
any image position is given by (Scheider et al. 1992, Section 
2.3) 



(8) 



where /i is the magnification. The critical curves in the image 
plane can be found by setting this equation equal to zero (/i — > 
oo), which yields two conditions: 
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The first condition produces a single caustic point at the ori- 
gin of the source plane, while the second condition produces 
two concentric caustic circles. In the image plane, the first 
condition is responsible for the outer of two critical curves, 
and the second condition for the inner plus a critical point at 
the origin. Figure Q] illustrates the transformation from the 
CIS potential to the PM+CIS potential in terms of the image 
positions and magnifications as a function (3. 



2.1.2. Point Mass and Broken Power Law Isothermal Sphere 
(PM+POW) 

We now derive a similar set of expressions for lensing by a 
point mass potential added to a circularly symmetric broken 
power law. In this case, the surface mass density is a more 
natural starting point and is given by ( Keeton 2003) 
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where 7 and 77 are inner and outer power law indices, respec- 
tively, 6b is the break radius between the power laws, I/, is 
the surface brightness at the break radius, a gives the sharp- 
ness of the bre ak, and T is the mass-to-light ratio. This yields 
(Keeton 2003) a lensing deflection angle, 
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where S cr is the critical surface density for lensing and 2^1 is 
the hypergeometric function. This model has introduced sev- 
eral new parameters that are not present in the cored isother- 
mal sphere. Later, we shall wish to explore the effects of these 
new parameters and compare the outcome back to the cored 
isothermal model. To facilitate that effort we now restrict the 
outer power index to its isothermal equivalent, r\ = 1, and rec- 
ognize that, under this condition, in the limit 8\, — > the sur- 
face mass density must match that for the singular isothermal 
sphere. Doing so, we have 
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and are able to write the deflection angle for the combined 
point mass and broken power law in terms of 6>sis: 
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2.1.3. Typical Values of 9pm, &c, $t» ct, and"/ 
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The empirical correlation in the local universe between the 
masses of SMBHs and the velocity dispersions of their host 
galaxy bulges can be exploited to determine a typical value 
of #pm/#sis- Observable lenses are dominated by elliptical 
lensing galaxies (Turner et al. 1984). For elliptical galaxies 
we can ignore the distinction between the velocity dispersion 
of the bulge and the galaxy as a whole. Therefore, we assign 
the SMBH its mass according to (Merritt & Ferrarese 2001): 



M BH = 1.48±0.24x 10 S M G 



200 km/sec 



4.7 



(17) 



TABLE 1 

Average values for the "Nuker Law" parameters 



Parameter 



Mean 



Median 



a 

7 

V 

9b (arcsec) 



2.35 
0.09 
1.28 
1.71 



2.03 
0.08 
1.33 
1.29 



NOTE. — These averages were calculated for the "cored" galax- 
ies listed in Faber et al. 1 1997, Table 2). Those authors use to 
indicate the outer power law index, whereas we have used rj. 



where a is equal to the (bulge) velocity dispersion and Mq is 
the solar mass. If we temporarily assume that this relation- 
ship holds for all galaxies, even beyond the local universe, we 
can generate characteristic curves for 6*pm/#sis as a function 
of <r, Dol, Dos, and Dis by substituting the above result into 
Equation|4] Figure |2] illustrates the behavior of 6*pm/0sis for 
different lens and source galaxy distances. For lenses located 
roughly halfway to the source, which because it is most effi- 
cient for lensing will dominate (Turner et al. 1984, Figure 5), 
there is not much variation in 6>pm/#sis f° r a fixed velocity 
dispersion. Galaxies with different velocity dispersions will 
produce different #pm/#sis ratios. 

There is to date no empirical understanding of any relation- 
ship between galaxy cores and other observables similar to 
that for SMBHs. The parameterizations we have chosen use 
angular measures for the core and break radii which, when 
mapped back to physical lengths, are dependent on the red- 
shifts of the lens and source galaxies, complicating their in- 
terpretations. Figure[3]gives the relationship between the an- 
gular and linear radii. (From here on we will simply refer to 
8q and 8b as the core radius regardless of the model under 
discussion.) 

For the PM+POW model, we have already restricted the 
outer power law index, rj = 1, to match the isothermal poten- 
tial, but left 7 and a unconstrained. Figure |4] illustrates the 
effects of varying these parameters. Drawing on the reported 
"Nuker Law" model fits, we can calculate average values for 
these parameters. In Table ^ we compute the mean and me- 
dian v alues for the 26 "cored" galaxies listed in Faber et al. 
( 1997, Table 2). 

Based on these findings, we will take 7 = 0.1 and a = 2.0 for 
all subsequent calculations. From Tabled it is also apparent 
that forcing 77 = 1, as we did above, was not a bad approxima- 
tion to these observational results. 

2. 1 .4. Shear and Ellipticity 

We constrain our analysis to circularly symmetric lensing 
potentials for this paper in order to better illustrate the de- 
pendence on fundamental parameters. As it is well known, 
however, gravitational lenses seem to rarely display such sim- 
ple symmetry and are often best modeled by elliptical mass 
density potentials or by including a gravitational shear in the 
local environment. In Appendix A we present a justification 
for neglecting the contribution of these higher order factors in 
our exploration of the fundamental behavior of these systems. 

2.2. Signatures of SMBH Lensing 

The addition of a SMBH to an otherwise typical galactic 
lensing potential produces three significant effects: 1) two- 
image lensing extends infinitely beyond the Einstein radius 



of the galaxy, 2) the maximum number of images increases 
from three to four, and 3) the area of the region in the source 
plane expected to produce more than two images is reduced 
compared to a cored or broken power law galactic potential 
with no SMBH. In general, for most source positions, (3, the 
lensing characteristics of the combined SMBH and galaxy act 
as if they were merely superimposed, and only conspire in the 
central regions of the source plane to cancel out the secondary 
images that either would have produced alone. 

Both the cored and broken power law potentials demon- 
strate these properties. In both models, the source plane is 
divided into three regions by two concentric circular caustics 
and a degenerate caustic point at the origin. In fact, as shown 
in Figure [2 with the PM+CIS model, for different values of 
the source position, either one or three secondary images may 
form opposite the position of the primary image. For the cases 
in which we are interested, the core radius and point mass 
have little effect on the properties of the primary image ("A"), 
but can dramatically change the properties of the three possi- 
ble secondary images ("B", "C", and "D"). 

Image B, where it exists, will be a saddle point in the light 
travel time and will have negative magnification. It will lie 
very close to the location of the secondary image that would 
form from a singular isothermal sphere potential alone, and 
therefore, provides little information about the properties of 
the core of the galaxy or the SMBH. Image C will be a maxi- 
mum of the light travel time and will have positive magnifica- 
tion. It will carry information about the degree of softening of 
the galaxy's potential, generally becoming brighter for larger 
cores and fainter for smaller cores. Image D will be another 
saddle point in the light travel time and again will have nega- 
tive magnification. For most source positions, it will track the 
position and magnification of the image that would form from 
an isolated point mass potential, thus providing information 
about the SMBH. Image D is present even without images B 
and C for source positions with sufficiently large (3. 

Keeton (2003) points out that including SMBHs in lensing 
calculations has the affect of "swallowing" faint core images. 
Figure ^ demonstrates that for sufficiently small impact pa- 
rameters, where core images would be faintest, the SMBH 
does indeed eliminate the core image that would have formed 
in its absence. Furthermore, Figures0]and|5]demonstrate that, 
for a given #pm/#sis, the s i ze of the region excluding core im- 
ages is dependent on the other properties (fie, or Ob, 7, and 
a) of the inner region of the host galaxy. Referring again to 
Figure 13 it is easily seen that for both models, reducing the 
core radius pushes the caustics out and closer together un- 
til ultimately they merge (this actually occurs before the core 
radius goes to zero, although the exact configuration is not 
shown in the plotted sequence). Therefore, smaller core radii 
decrease the likelihood of observing four-image lenses as can 
be seen in Figure [6] where we have plotted the area between 
the caustics as a function of core radius for multiple curves of 
constant 6*pm/#sis- This plot also demonstrates that increas- 
ing #pm/#sis also acts to reduce the region of four images. 
This degeneracy accounts for the weake ning o f constraints on 
inner mass-profiles discussed by Rusin & Ma (2001). 

3. NUMERICAL ANALYSIS OF THE MODELS 

3.1. Lensing cross-sections 

As we have seen, the PM+CIS and PM+POW models pro- 
duce lenses with a minimum of two images and a maximum of 
four images, depending on the source position. We would like 
now to extend the analysis and gain insight into the likelihood 



of observing these images, particularly image C and image D. 
To do so, we must introduce a parameter which relates to the 
observational ability to detect an image. Let us assume, for 
simplicity, that an image can be detected if it produces a flux 
that is above a specified fraction of the flux produced by the 
primary image: 



Fb,c,d _ Hb,c,d > _J_ 
F A (j, A ~ R : 



(18) 



where R is the maximum detectable dynamic range between 
images and is the parameter we will use in the following cal- 
culations. The cross-section for detection is then defined as 
the area in the source plane that produces images which meet 
this requirement, and is in general characterized by two ex- 
treme source positions, (3 m { n and /3 max , such that the total an- 
gular area, A, is given by 

A = n((3 2 mRX -(3l in ). (19) 

In this expression, A, /?„„„, and j3 max are implicitly dependent 
on the core radius, the SMBH Einstein radius, the galaxy Ein- 
stein radius, and the maximum dynamic range. 

We are now in a position to calculate the characteristic lens- 
ing cross-sections for the PM+CIS and PM+POW models. 
We consider three cross-sections: 1) for detecting two im- 
ages (A and B, or A and D), 2) for detecting three images 
(A, B, and C), and 3) for detecting four images (A, B, C, and 
D). The cross-section for observing three images is, in fact, 
the sum of two of the annuli discussed above, one along the 
inner caustic and the other along the outer caustic. Due to 
the nature of the lensing potential equations, numerical meth- 
ods were employed to determine /3 m ; n and /3 max for a range of 
model parameters to a precision of 0.002 6*sis- The results are 
discussed below. 

3.1.1. Two Images 

It is clear from Figure [5] that combining a point mass and 
isothermal sphere produces the possibility for a very large 
two-image cross-section. In fact, without a maximum dy- 
namic range constraint the cross-section would be theoret- 
ically unbounded, and even for finite dynamic range con- 
straints, /3 max may be well beyond the standard singular 
isothermal sphere limit of /3 max = 6*sis- To see more quanti- 
tatively how dynamic range affects the cross-section, Figure 
0plots At as a function of core radius for several values of R, 
ranging from 1 to 10 8 , with 6*pm/#sis fixed at 0.03. There are 
two striking features in this plot. First, there is a degeneracy 
of cross-section curves, around A2/(tt01 is ) ~ 1 (the same fea- 
ture is visible in Figure^). This is caused by the disparity in 
magnification between image B and the image D tail, both of 
which can contribute to the two-image cross-section. The dis- 
parity reduces as the core radius is reduced and eventually, for 
cores less than about 6*c/#sis — 10~ 2 8 , disappears because the 
source plane is no longer divided into three regions, but only 
one that always produces two images (see Figure [5] particu- 
larly the bottom row). This situation eliminates the distinc- 
tion between our labels "B" and "D" for images, although, 
for (3 < 6*sis the secondary image has properties similar to 
those of a typical image B and, for (3 > #sis, the magnification 
drops toward that of a typical image D. This smooth decrease 
in magnification near 6*sis for lenses with small core radii is 
responsible for the second noticeable feature: cross-sections 
increase rapidly with dynamic range for small core radii since 
there is no step function in magnification to overcome as is 
the case for larger core radii. 



The SMBH has an effect on the cross-section, as well. Fig- 
ure [8] specifically addresses this dependence, and from it we 
see that the two-image cross-section increases with #pm/#sis- 
This is due to the image D tail increasing in magnification for 
larger SMBH mass scales. The effect is most dramatic for 
small core radii. 

3.1.2. Three Images 

For the three-image cross-section, it is also useful to begin 
by examining the limit of no dynamic range constraints. In 
this case, the region in the source plane that produces three 
observable images is defined completely by the two concen- 
tric caustics as discussed at the end of Section 2 and plotted in 
Figure[6] This represents an upper bound for the three-image 
cross-section which is plotted in Figures [8] and Large dy- 
namic ranges will saturate to this limit. As the dynamic range 
is lowered, the cross-section reduces with the most rapid de- 
cline at small core radii. This can be understood from exam- 
ination of Figure |5J once again noting that the magnification 
of the image C generally decreases as core radius decreases. 
It is no coincidence that the range of core radii that produce a 
significant three-image cross-section corresponds to the range 
where the two-image cross-section is highly degenerate since 
both trace the portion of parameter space that is dominated 
by the region in the source plane between the two concentric 
caustics. 

Unlike the two-image cross-section, increasing the SMBH 
mass scale decreases the three-image cross-section. This is 
another manifestation of the image "swallowing" previously 
described in Section|2] 

3.1.3. Four Images 

The four-image cross-section is qualitatively similar to the 
three-image cross-section. However, since the overall mag- 
nification of image D is considerable less, the cross-section 
curves are correspondingly reduced as is clear in Figures [8] 
and|7] 

3.2. Lensing Statistics 

In this section we derive the optical depth to lensing 
for a given source using the calculated cross-sections from 
above and following the standard method (Turner et al. 1984; 
Scheideret al. 1992, Chapters 11 and 12). If we assume that 
the lensing objects have a local comoving number density that 
can be written as «/,(cr,z), then the differential optical depth of 
a lensing event occurring at a redshift zl for a source at zs is 



dr ■■ 



■■n L (a,ZL)(l+ZL) 3 A(a,ZL,Zs)Do L (zL)dacdt. 



(20) 



Here we constrain #pm/#sis according to Equation 1171 and 
have explicitly written the cross-section, A, as a function of 
velocity dispersion, lens redshift, and source redshift to high- 
light its dependence on these parameters. It is, of course, still 
dependent on the choice of core radius and dynamic range, as 
well. Dol is again the angular diameter distance between the 
observer and lens, an d, in a homogeneous, isotropic universe 
we have llHo ggI1999l) 

A " [n M (l+ZL) 3 +n k (l+Z L ) 2 + n A ] l/2 dZL, (21) 



cdt ■ 



l+z L 



where c is the speed of light, t is the lookback time, and 
Rq, VIm, ^a> and $7* are the standard cosmological parame- 
ters with ft k = 1-£Im-Qa- Conceptually, all that remains 
is to select the form of the density distribution function. It 



TABLE 2 
Luminosity function parameters for spiral, SO, and 

ELLIPTICAL GALAXY TYPES 



Parameter 


Spiral 




SO 


Elliptical 


a 


-1.16 




-0.54 


-0.54 


n„ 


1.46 X 1(T 2 


0.61 x 1(T 2 


0.39 X 1(T 2 


7 


2.6 




4.0 


4.0 


(T* 


144 




206 


225 
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is common in gravitational lens studies to utilize either the 
Press-Schechter formalism (Press & Schechter 1974) or the 
Schechter luminosity function (Schechter 1976) to generate 
this distribution. We have chosen to use the Schechter lumi- 
nosity function combined with the Faber-Jackson relation and 
have also assumed no number density evolution with redshift 
to cast the density distribution in terms only of velocity dis- 
persion giving: 



«l(ct) = 7«* — 
\<J* 



7(a+l)-l 



exp 



(22) 



We use the same choices of parameters a, 7, and er*, and 
abundances , n», o f spiral, SO, and elliptical galaxytypes as 
Ofek et al. (2003) and have reprinted them in Table [2] 

The redshift distribution function of extragalactic flat- 
spectrum radio sources is localized around z s — 1 and has a 
mean value of (z s ) = 1.1 (Munoz et al. 2003). We will assume 
that treating all sources as if they were at the mean redshift 
is sufficient and take z s = 1.1, We are now prepared to cal- 
culated the optical depth to lensing for each cross-section of 
interest by integrating over velocity dispersion and redshift. 
Again these calculations must be carried out numerically be- 
cause the cross-sections themselves are not analytical. 

The results of these calculations are presented in Figure [9] 
For two-image lensing, we see that the optical depth to lens- 
ing, t, parallels the cross-section curves for this case. For 
high maximum dynamic range, the probability is quite high, 
and as the dynamic range is lowered to more typical values 
of R ~ 1000, the probability approaches the expected value 
for lensing by a population of singular isothermal spheres of 
t ~ 0.001 (Browne et al. 2003), as it must since this result has 
been observationally justified. 

The optical depth to lensing for three- and four-image cases 
are also similar to their cross-sections, with the four image 
case suffering a more rapid decline in optical depth as R is 
decreased. Both probabilities exhibit significant dependence 
on the assumed core radius of the lensing population, indicat- 
ing that the addition of a point mass to a softened isothermal 
lens mode l, if restricted to the relationship discussed in Sec- 
tion l2.1.3l would not hinder using lens statistics to learn about 
galaxy cores. 

4. EMPIRICAL IMPLICATIONS OF THE MODELS 

4. 1 . Relevant Facilities and Instruments 

We have seen that if an intervening lens galaxy contains a 
SMBH, then an additional faint image is predicted due to the 
black hole potential which is not present in unmodified mod- 
els. In fact, even in cases when the source position and the 
lens are too far apart to create multiple images in the standard 



way, the SMBH will still produce a faint image of the back- 
ground object near the center of the foreground galaxy. In 
this section we discuss the applicability of optical, infrared, 
and radio facilities to observing these images. 

4. 1 . 1 . Optical and Infrared 

Optical observations of these highly demagnified images 
would be difficult. In order to obtain a significant cross- 
section for detecting image D, we saw in Section[3land Figure 
01 that demagnification from the primary image by 10 , or 15 
magnitudes, must be observable. In the Hamburg-ESO Sur- 
vey (Wisotzki et al. 2000) and the Sloan Digital Sky Survey 
(SDSS), the brightest quasars at redshift greater than 0.5 have 
magnitudes of 15. The demagnification due to the lensing 
alone would reduce this to 30th magnitude, and this would 
be compounded with the extinction that the image would suf- 
fer passing through the inner regions of the lensing galaxy. 
Such faint observations are far beyond the capabilities of cur- 
rent optical telescopes. To observe the lensing of even the 
brightest quasars would require 30-meter class or larger opti- 
cal telescopes. 

Infrared observations would not suffer the severe extinction 
of the image passing through the foreground galaxy, however, 
it would still be problematic to separate the lens image from 
foreground galaxy emission. The use of spectral lines to dis- 
tinguish two redshift components may be able mitigate this 
problem. 

4.1.2. Radio 

Observations with radio telescopes, while still difficult, are 
more feasible. The brightest sources above a reds hift of 
z = 0.5 in the NRAO-VL A Sky Survey (Condon et al. 1998) 
and the FIRST sky survey (Becker et al. 1995) have flux den- 
sities greater than 20 Jy. The flux density of image D from 
such a source would be on the order of 20 /iJy, and at ra- 
dio wavelengths any absorption through the lensing galaxy 
should be negligible. Long integrations with the VLA and 
VLBI arrays can already achieve sensitivities of 10 /iJy per 
beam (Garrett et al. 2003) or better, sufficient to observe faint 
images from the very brightest radio sources. Future improve- 
ments in radio telescope instrumentation, however, will dra- 
matically increase the predicted number of observable D im- 
ages. Four projects in particular are relevant to making sensi- 
tive observations of lenses and are summarized below. 

VLA Expansion Project — The VLA Expansion Project is al- 
ready underway to supply the VLA with modern electronics 
and also plans to add up to eight additional antennas to the 
array at distances out to 350 km from the core facility. The 
resulting Expanded VLA (ELVA), scheduled for completion 
by 2012, is intended to improve performance in important ar- 
eas by a factor of 10. This is expected to extend the detection 
threshold of faint sources to 1 /iJy or less. With this sensitiv- 
ity, the EVLA should also be able to detect the C images for 
existing radio lenses. 

Square Kilometer Array — The Square Kilometer Array (SKA) 
is a proposed next generation radio telescope that offers the 
most promise for observing lenses with faint images. The cur- 
rent design goals for the SKA are to achieve a dynamic range 
approaching 10 8 and a detection threshold that could reach to 
the 10 nJy level. 



Low Frequency Array — A low frequency radio telescope ar- 
ray has been proposed that would operate at low frequencies 
between 30 and 240 MHz. The design is intended to achieve 
a baseline dynamic range of 10 5 , extending to 10 8 in some 
cases, and a detection threshold of less than 1 /iJy. However, 
due to the low frequencies used, the expected best-case angu- 
lar resolution of about 1 arcsec (at 240 MHz) of the instrument 
would prevent some lens configurations from being resolved. 

Atacama Large Millimeter Array — The Atacama Large Mil- 
limeter Array (ALMA) is a millimeter wavelength telescope 
under design by Europe and North America. It will be located 
in Llano de Chagnantor, Chile, and consist of at least 64, 12- 
meter antennas. The intended resolution is 10 mas wit h a sen- 
sitivit y approximately 20 times greater than the VLA. (Tofani 
2003) 

4.2. Lensing Expectations 

Despite the limitations of current technology, it is not too 
early to consider candidate SMBH lensing targets. The first 
step is to estimate the predicted total number of lenses that 
should be observable. We can use the optical depth calcu- 
lations from the previous section along with the known flux 
distribution of radio sources to accomplish this. To estimate 
the number of lenses an instrument with a limiting sensitiv- 
ity, S, should be able to detect, we integrate the product of the 
differential source flux distribution, dN/dF, and the optical 
depth, t(R,6c). The appropriate dynamic range, R, is given 
by the ratio between the sensitivity of the instrument and the 
flux of the source (R = F/S): 



N(S,8c)-- 



(F/S, 



dN_ 
d~F 



dF. 
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We are neglecting the dependence on source redshift of both 
t and dN/dF. Again, we will assume that treating all sources 
as if they were at the mean redshift, z, = 1.1, is sufficient. We 
still cannot perform this integral fully because we have only 
evaluated r for specific values of R. However, we can achieve 
an approximate evaluation by converting the integral to a sum 
over flux bins centered on equivalent R values: 
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The differential number of radio sources as a function of flux 
in the range 14.5 < F < 1500 ^zJy is (Windhorst et al. 1993): 



dN 

dF^~ 



= (-4.6±0.7)F 



-2.3±0.02 
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Jy-'Sr" 
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We will use this relation for all flux densities, although it is 
likely to overestimate the the number of faint sources (F < 1 
nJy) and underestimate th e number of bright sources (F > 10 
mJy) as discussed by Windhorst et al. ( 1993). Finally, in or- 
der to perform the sum, we must choose a core radius and an 
image of interest. Two diagnostic core radii regimes are evi- 
dent from Section|3] We begin by considering: 1) the limit of 
vanishing cores, and 2) the core radius that produces the max- 
imum likelihood of observing three- and four-image lenses. 
The first regime is approached by setting #c/#sis = 10~ 4 and 
performing the sums described above. Using the PM+CIS 
model we calculate, for each image, the number of observ- 
able lenses per steradian. 

Table|3]illustrates that no lenses with more than two images 
are expected. This is readily explained by noting in Figure 



TABLE 3 

Predicted number of gravitational lenses per 

steradian using the empirical flux distribution of 

KNOWN RADIO SOURCES AND 6>c/0SIS = Iff" 



TABLE 5 
Number of sources in the FIRST radio catalog 



n-4 



NOTE. — This table indicates the predicted number of gravi- 
tational lenses per steradian for a search with sensitivity, 5, using 
the optical depth calculations for the PM+CIS potential from the 
previous section and the empirical flux distribution of known radio 
sources. The choice of core radius is $c/$sis = 10 4 to approxi- 
mate the limiting case of no galaxy cores. Nj ota i is the number of 
two-image lenses predicted, N 3 is the number of lenses in the sub- 
set where a third, core image is detectable, and N4 is the number 
of lenses in the subset with a fourth, SMBH image. The last col- 
umn, labelled No uts jj,,, is the number of two-image lenses where the 
fainter image is due directly to the SMBH. 



TABLE 4 

Predicted number of gravitational lenses per 

steradian using the empirical flux distribution of 

known radio sources and 0c/#sis = 1(t l3 



S (ply) 


Example 


Nrotal 


N 3 


^4 


NoutsMe 


1000 


FIRST Survey 


4 











100 


Typical VLA 


73 


1 








10 




1,400 


18 








1 


EVLA 


29,000 


360 


1 


1 


0.1 




580,000 


7,100 


13 


16 


0.01 


SKA 


12,000,000 


140,000 


260 


330 



NOTE. — This table indicates the predicted number of gravita- 
tional lenses per steradian for a search with sensitivity, S. The choice 
of core radius is #c/#srs = 10 3 , which is the value of 6q that pro- 
duces the largest optical depths for core and SMBH images. 



[9] that the optical depths for three- and four-image lenses at 
#c/#srs = 10~ 4 are negligible for all flux ratios. The last col- 
umn in the table, labelled "Noutside", indicates the number of 
two-image lenses where the secondary image would not be 
present without the SMBH and, therefore, may still be able to 
provide useful information about the SMBH. 

The second core radius regime is that of maximum like- 
lihood to produce three- and four-image lenses. For the 
PM+CIS model, this occurs near 6> c /6>sis = 10" 1 3 . The cal- 
culated values for this core radius are listed in Table|4] 

These values indicate that if typical core radii for galaxies 
are in this range, we would expect to find, with an SKA-like 
instrument, no more than about 260 lenses per steradian with 
a fourth, SMBH image out of about 12 million. On the other 
hand, 140,000 lenses per steradian with a third image bright 
enough to detect would be expected. 

4.2.1. Limiting Candidate Targets 

Despite the tremendous number of lenses capable of pro- 
viding information about an SMBH that an extremely sensi- 
tive search is predicted to detect, it would be operationally 
infeasible due to the required long integrations times to con- 
duct such a campaign without first reducing the number of 



Bin 



Flux (mjy) 



N, 



bin 



S fjuJy) 


Example 


Nrotal 


N 3 


N4 


Noutside 


8 
7 


1000 


FIRST Survey 


4 











6 


100 


Typical VLA 


74 











5 


10 




1,500 








2 


4 


1 


EVLA 


30,000 








34 


3 
2 
1 


0.1 




590,000 








670 


0.01 


SKA 


12,000,000 








13,000 



10" 
10 30 

10 2.5 
10 2.0 

10 L5 
10'° 

10 05 

1 



27 

227 

1,444 

7,418 

26,914 

75,756 

181,766 

517,565 



NOTE. — The number of sources in the FIRST radio catalog are 
divided into 8 bins by flux. The flux values listed in the table repre- 
sent the floor of each bin. The number of sources in each bin is used 
in combination with optical depth calculations for the PM+CIS lens- 
ing potential to calculate the predicted number of lenses contained 
in the catalog (Tablesl6landl7V. 



TABLE 6 
Predicted number of gravitational lenses in the 

FIRST CATALOG USING C /0SIS = 10" 4 



S ( M Jy) 


Example 


N Total 


N 3 


N 4 


^Outside 


1000 


FIRST Survey 


100 











100 


Typical VLA 


470 








4 


10 


- 


650 








40 


1 


EVLA 


790 








180 


0.1 


- 


1100 








460 


0.01 


SKA 


1700 








1100 



NOTE. — This table indicates the predicted number of gravita- 
tional lenses for a search with sensitivity, S, in the FIRST catalog of 
radio sources. The choice of core radius is 9r/9sis = 10~ 4 . 



candidate targets. 

If only sources in a radio catalog are under consideration as 
possible targets, the sensitivity limit of the catalog will restrict 
the number of two-image lenses and increase the overall ratio 
of three- and four-image lenses to the total number of lenses. 

There are over 800,000 sources in the FIRST radio catalog 
(White et alJl!99 7) which has a detection limit of 1 mjy and 
covers nearly 10,000 square degrees of the sky, mostly in the 
northern hemisphere. For this analysis, we will assume as we 
did earlier, that all the sources in the catalog are at redshift, 
z= 1.1, and follow the same procedure as above. We divide 
the catalog into eight flux bins, shown in Table[5] 

and perform the sum for the two diagnostic core radii 
regimes. 

The first regime, with 9 C /0 SIS = 10" 4 and using the PM+CIS 
model, predicts the number of observable lenses for the area 
of the sky covered by the FIRST survey given in Table|6] 

Immediately it is clear that, compared to Table [3] the num- 
ber of total lenses detectable by sensitive instruments has been 
drastically reduced. This is because the FIRST survey con- 
tains no sources with flux below F = 1 mjy and, therefore, the 
extra sensitivity has no accessible faint sources on which to 
operate. 

The second core radius regime of interest is where Oc/Osis = 
10~ 1,3 . The calculated values for this core radius are in Table 

m 

Although significantly fewer lenses are likely to be con- 
tained in the catalog, the relative number with an observable 



TABLE 7 

Predicted number of gravitational lenses in the 

first catalog of using 6» c /6»sis = 10 _1 3 



S (AJy) 


Example 


Nrotal 


N 3 


N 4 


Noutside 


1000 


FIRST Survey 


90 


2 








100 


Typical VLA 


410 


20 








10 


- 


410 


130 


1 





1 


EVLA 


410 


320 


12 


1 


0.1 


- 


430 


320 


70 


7 


0.01 


SKA 


580 


320 


230 


86 



NOTE. — This table indicates the predicted number of gravita- 
tional lenses for a search with sensitivity, 5, in the FIRST catalog of 
radio sources. The choice of core radius is #c/#srs = 10 -1 ' 3 . 



core image is much higher (between 50% and 80% for EVLA 
or better levels). In addition, most of the two-image lenses 
should already be detectable with the VLA. It is plausible, un- 
der these circumstances, that an appropriate search guided by 
the FIRST catalog would be able to detect and characterize, 
today, a population of two-image lenses that will yield a sig- 
nificant number of third and fourth images when re-observed 
with future technology. 

The above calculations were performed for two specific 
core radii of particular interest, but they can be repeated for 
all the core radii we considered in Section [5] This has been 
done and is plotted in Figure^] for S=l /zJy. The ratios of 
N-i/Nj t a i and Na/NjoioI are plotted as well. For sensitivities 
better than 10 fiiy, the A^/A^om/ increases rapidly with core 
radius in the region ^c/^sis > 10~ 26 and approaches unity for 
#c/0sis = lO" 06 and S=10 /Jy. Similarly, N4/N To tai increases 
in the same region, although it reaches a maximum of only 
approximately Na/NtoiuI — 0.4. 

Searching for 500 lenses out of a catalog of over 800,000 
sources is still a daunting task. We suggest that another object 
catalog may help reduce the list of candidates even more. The 
Sloan Digital Sky Survey (SDSS) covers 3324 square degrees 
as of Data Release 2, overlapping the FIRST survey, and will 
eventually cover nearly 10,000 square degrees. It has identi- 
fied over 88 million objects at optical wavelengths, many of 
which are extragalactic sourc es with est imated redshifts be- 
tween < z < 0.5 (Stoughto n et al.ll2002l) . This population 
of extragalactic objects makes an ideal foreground of poten- 
tial lens galaxies for background radio sources in the FIRST 
catalog. Ilvezic et all (120021) have shown that the error in posi- 
tion between the two surveys is ~ 0.4 arcsec per coordinate, 
and the SDSS has included in its database a special flag to in- 
dicate objects that are spatially collocated with FIRST radio 
sources.. If pairs of neighboring objects within approximately 
1 arcsec of each other could be found in the SDSS database 
where one is marked with the FIRST flag, then this would 
represent a sample of targets with an elevated likelihood of 
containing radio lenses. If it could be further reduced to pairs 
where the radio source is known to be behind the neighboring 
object, the likelihood would be higher, yet. Fortuitously, the 
SDSS aims to acquire 1 million spectra of galaxies, preferen- 
tially targeting identified optical counterparts of FIRST radio 
sources. Thus, redshift information may be available for one 
or both of the objects in a pair directly from the database. 

A complimentary approach to finding new lens candidates 
has been described by Bolton et al. (2004). They have been 
able to identify multiple emission lines at two redshifts in a 



number of SDSS spectra, indicating a foreground and back- 
ground galaxy within the 3 arcsec field of view, and are using 
these objects in an optical lens search. 

4.2.2. Galaxy Cores 

It is clear from the preceding estimates that the size of 
galaxy cores greatly affects the likelihood of observing core 
and SMBH images. It is beneficial, therefore, to consider 
what is known about galaxy cores to date. The structure of 
the central regions of galaxies has been studied extensively 
through surface brightness measurements and a number of at- 
tempts have already been made to deduce the properties of 
galaxy cores from gravitational lensing studies. 

Fits to the "Nuker Law" surface brightness models, which 
we used earlier as the basis of our second lens model, have 
given the result that galaxies appear to be of two types, those 
that are "cored" and those that are not. The cored galaxies are 
those that are luminous, slowly rotating ellipticals with boxy 
isophotes; the non-cored, power law galaxies are faint, rapidly 
rotating galaxies with disky isophotes. There is a fairly tight 
correlation between break radius and luminosity. Large el- 
liptical galaxies dominate the lensing cross-section, implying 
that the galaxies that are important for SMBH lensing are also 
those that have cores. In fact, it has been hypothesized that in 
the hierarchical merger process involved in galaxy formation, 
SMBH binary systems will form and spiral toward the center 
of the new system, cooling and possibly ejecting starts near 
the center. The SMBHs may, therefore, even be responsible 
for the existence of the cores. 

The presence of a core image (image C) in grav- 
itational lenses, however, is exceedingly rare, with 
candi date images observed only in PMN J1632-003 3 
dWinnetaU 12001 . MG 1131+0456 (iChen & Hewitt! 11991 . 
and APM 08279+5255 lllbata et al.lll999l) . For thelatter two 
systems, it remains uncertain whether third component is an 
image or emission from the lens galaxy. If we do take the 
third components in MG 1131+0456 and APM 8279+5255 
as lens images, best fit values of 6WflsT,s p ^ 0.2 (Ch en et all 
119951) and 6» c /6>sis - 0.7 ( lEgami et alJ20 00). respectively, are 
calculated. 

The only system with convincing evidence of a c ore im- 
age detection is PMN J1632-0033 (Winn et aT] l2004l) . Radio 
spectra of all three image components have been acquired and 
the flux ratios Fb/Fa and FqJFa are nearly constant across the 
high frequency measurements indicating the components are 
likely from the same source. The only discrepancy is with the 
lowest frequency measurement at 1.7 GHz, but as Winn et al. 
(2004) argue, this is where absorption from passing through 
the dense part of the lens galaxy would be strongest. They 
have, furthermore, been able to place constraints on sev- 
eral quantities of interest (Winn et al. 2003). Using a broken 
power law model, the relation between break radius and inner 
power law slope is constrained; and using a cored isother- 
mal mass-density model, a core radius of 8c — 5.5 mas or 
r c ~ 31/T 1 pc, is derived. Compared to the estimated Einstein 
radius of the lens galaxy of #sis — 0.74 arcsec (Winnet al. 
2002), the core radius corresponds to #c/#sis = 0.0075, al- 
though this is not directly comparable to our core radius due 
to the differing model. 

Keeton (2003) discusses in detail the implications of the be- 
havior of the mass-density profile at the core to the detectabil- 
ity of a third image, and concludes that the paucity of third 
images in lens systems is co nsistent with observed surface 
density profiles. Rusin & Ma (2001) have studied six double 



lens systems in which the third image is not detected. Using 
the known detection limits for each of these sources, they em- 
ploy a singular mass profile with variable slope to determine 
how steep the inner density profile must be to explain the lack 
of this image. They find that for these six sources the inner 
region must be quite steep. However, it is not clear that the 
constraints on the inner slope are relevant in constraining the 
possible flattening of the mass distribution that might occur 
right at the SMBH. 

The above attempts have tended to be specific case studies 
of individual lenses and therefore it is difficult to get a sense 
for the overall probability of any particular core radius. The 
results of our modelling are in a position then, to provide the 
backdrop for interpreting these results. For example the re- 
sults of a systematic lens survey with 1 mJy sensitivity could 
be compared to the ratio of N-s/Njotai m Figure ^3 to learn 
about typical core radii values. 

4.3. Structure in Image D 

To this point we have considered only the overall magnifi- 
cation and position of image D. If these images are detected, 
however, any observed structure in them will be capable of 
providing additional important information about the lens sys- 
tem. The properties of observed structure in a lens image 
should be distinct from the properties of any radio emission 
from the foreground galaxy, and thus, resolved structure in 
faint emission at the center of a lens could be used to reduce 
confusion as to whether an image or foreground emission had 
observed. Additionally, structure in image D is capable of in- 
dicating the presence of binary SMBHs in the lensing galaxy. 
Both of these are discussed in detail below, but first we will 
examine the ability to resolve structure in these images. 

4.3.1. Detecting Structure in Image D 

The flux density of image D, as we have seen, is demagni- 
fied typically by a factor of at least 10 6 from that of the orig- 
inal source. Since surface brightness is conserved by lensing, 
this corresponds to a reduction in the surface area of the image 
compared to the original source. Therefore, we expect each 
of the two spatial dimensions on the sky to be reduced by a 
factor of order VT0<> = 10 3 . 

VLBI measurements achieve angular resolutions approach- 
ing 0.1 mas. Thus a radio source with 1 arcsec scale struc- 
ture would produce an image D with ~ 1 mas scale structure, 
which should easily be distinguished by the VLBI. 

4.3.2. Foreground Galaxy Emission 

Observing a faint signal near the center of a lens galaxy 
does not provide sufficient evidence to prove that image D has 
been detected. It is possible that the signal may be explained 
by a lens model without a SMBH or that the emission may be 
from the lens galaxy itself. 

In the former case, an analysis for plausible lens models of 
the probability of producing a faint image close to the cen- 
ter of the lens galaxy should typically be sufficient to remove 
confusion. For example, a standard singular isothermal sphere 
model produces a faint secondary image in rare cases where 
the source impact parameter is just inside the Einstein radius. 
This image, in principle, could be confused with the image 
D tail of the models we have discussed above. However, the 
probability of such a configuration in the SIS model is ex- 
tremely low compared to that for the PM+CIS and PM+POW 
models. 



In the latter case, structure in the signal could be used to rule 
out the possibility that the detected signal is emission from the 
foreground galaxy. Gravitational lenses produce a character- 
istic arcing in images that are resolved. In addition, structure 
in secondary images should be correlated with structure in the 
primary image. 

As an illustrative example of the lensing that we expect to 
occur due to a SMBH within a galaxy, we consider a radio-jet 
QSO that is known to have nearby foreground galaxies. The 
radio and optical observations of 3C9 are discussed in detail 
by Kronberg et al. (1991). This source lies at a redshift of 
Zs = 2.012, has total flux of approximately 190 mJy, and its 
radio jets are approximately 10 arcsec across. A radio map 
of the source is shown in Figure ^2 which was constructed 
from data provided by P. Kronberg. The especially interesting 
feature of 3C9 is that there are two foreground galaxies near 
the jets. One (galaxy 1) lies approximately 10 arcsec from 
the jets, has an R-band magnitude of 19.5, and has a mea- 
sured redshift of zl = 0.254. The other (galaxy 2) lies only 
about 3 arcsec from the jets, has a R-band magnitude of 23.0, 
and has a less well determined redshift. We will assume that 
the two galaxies lie at the same redshift (zl = 0.254). Using 
a simple ray tracing algorithm, we simulated the lensing that 
would occur from a SMBH in the center of each galaxy. Since 
the foreground galaxies are far from the radio source, and in 
such cases the properties of image D are determined almost 
entirely by the SMBH without influence from the host galaxy 
potential, we are able to work in the simplifying approxima- 
tion of treating the SMBH as an isolated point mass. Figure 
1 121 shows the image that would result from a 10 9 M Q SMBH 
in the center of galaxy 1 . The total flux from this image is ap- 
proximately 0.03 nJy. A SMBH with the same mass in galaxy 

2 would produce a much brighter image with total flux about 

3 nJy. Figure[T3lshows this case. 

These faint images of 3C9 would not be detectable, even 
with the anticipated sensitivity of the SKA, and are presented 
merely as examples of how structure could be identified. 3C9, 
therefore, is not a realistic candidate for future SMBH lensing 
studies. 

4.3.3. Mergers and Binary SMBHs 

An interesting question arising from the hierarchical for- 
mation of SMBHs is how long the black holes of two merg- 
ing galaxies take to coalesce. Indeed, if the time scale is 
long enough we would expect some galaxies to contain bi- 
nary SMBHs. Yu (2002) has carried out theoretical studies 
that determine how these merger time scales vary according 
to the type of host galaxy, its velocity dispersion, and the 
mass ratios of the SMBHs involved. The "surviving" binary 
systems (those whose merger time is larger than the Hubble 
time) are found to occur predominantly in galaxies that are 
not highly flattened and that have dispersion velocities greater 
than 90 km/s. In addition, SMBH pairs with equal masses are 
more likely to survive. 

The separations of most of the remaining binary systems 
range from 0.01 to 10 pc. From Equation [4] we calculate 
that for SMBHs with masses between 10 6 and 10 9 M Q , typ- 
ical values of #pm correspond to physical distances between 
10 and 400 pc. The binary separations, therefore, may be 
significant fractions of the SMBH Einstein radii and binaries 
may be distinguishable in some lenses. Binary SMBHs with 
the widest anticipated separations may be sufficient to alter 
not only properties of image D, but also image C (and there- 
fore the lensing statistics of these images), because for small 



10 



source impact parameters image C is predicted to form typi- 
cally within about 10 pc of the center of the lens (see Figures 
[Jand[3}. For large source impact parameters, however, image 
C does not exist and only image D will be affected. 

Figure fT4l demonstrates the results of considering two 5 x 
10 8 M© SMBHs in galaxy 2 of the 3C9 example consid- 
ered above. We calculate for these SMBHs that the Einstein 
radii correspond to 250 pc (using Q,m = 0.3, il\ = 0.7, and 
Hq = 70 km/s). It is evident that, even for separations of a 
few parsecs, noticeable effects are visible in this scenario. A 
survey of a large number of these faint lenses would provide 
insight not only into the distribution of the masses of SMBHs, 
but also whether the SMBHs occur as individual or binary 
systems and their correlation to the morphological type of the 
host galaxy. 

It is also evident from theoretical simulations that, either 
through mergers or initial conditions, individual SMBHs may 
not be located precisely at the center of their host galaxies. 
If a SMBH has formed at a distance greater than approxi- 
mately 100 pc (or more, depending on its mass) from the 
center of its host galaxy, Ullio et al. (2001, Figure 3) have 
predicted that dynamical friction is insufficient to produce 
spiral-in timescales shorter than the Hubble time. Addition- 
ally, Brownian motion may cause a SMBH or binary pair to 
wander around the nucleus of a galaxy within a radius of order 
0.01 pc (Merritt 2001). Observational evidence, on the other 
hand, suggests that in nearb y galaxies SMBHs are within 
0.1 pc of galaxy centers (Milosavlievic & Merritt 2001, Sec- 
tion 1). As was the case for binary systems, non-centrally lo- 
cated SMBHs may affect the magnifications curves for small 
source impact parameters, and thus the lensing statistics for 
images C and D. However, they are not likely to signifi- 
cantly alter the properties of systems with large source im- 
pact parameters since, in this limit, the properties of image 
D approa ch tho se of an image from an isolated SMBH (see 
Appendix IA.3I for a justification). Therefore, an off-center 
SMBH should not complicate detection or identification of a 
SMBH image in a two-image system with large source impact 
parameter. 

5. CONCLUSION 

Observational measurements of the relationship between 
SMBHs and the properties of their host galaxies are an impor- 
tant method for probing theoretical hierarchical growth mod- 
els. Gravitational lensing is a unique mechanism for acquiring 
the this information in systems at cosmologically significant 
redshifts. 

We have reviewed the calculations required to transform 
historically popular galactic lens models in order to account 
for the presence of SMBHs. The modified models, PM+CIS 
and PM+POW, exhibit an inner disk in the source plane where 
only two bright images are produced, surrounded by an anu- 
lus where two additional faint images are added near the cen- 
ter of the lens for a total of four images. For all source po- 
sitions outside this region, the models predict the presence of 
one primary image and one highly demagnified image deep in 
the center of the lens which traces very nearly the secondary 
image of an isolated SMBH. The relative cross-sections of 
these three regions are very sensitive to galaxy core sizes 
and SMBH masses. We used these characteristic signatures 
to demonstrate the origi n of the SMBH lensing effec ts re- 
porte d previously in lRusin & MalJ2001l) . lKeetonU2003h . and 
Maoetal1 d2(K)lli . 

Determination of lensing cross-sections over a range of pa- 



rameter space and subsequent optical depth calculations re- 
vealed that, even for the most favorable lens configurations, 
observations of the highly demagnified fourth images will re- 
main elusive until the next generation of radio telescopes are 
in use. Two scenarios, however, offer more immediate results. 
Since core images are detectable now, if galaxies tend to have 
large cores, deep VLBI observations (or EVLA) should make 
accessible a sufficiently large population of core images to 
provide meaningful statistics. On the other hand, if galaxies 
tend to have very small cores, then two-image lenses where 
the secondary image is the highly demagnified SMBH im- 
age already may be observable. For any scenario, there is 
strong evidence suggesting that a targeted search today for 
new lenses, utilizing sky surveys such as FIRST to limit the 
number of candidate objects, should yield a precursory sam- 
ple of two-image lenses that will, with improved technology, 
provide a significant number of third and fourth images. And 
finally, if these images are detected they may even reveal bi- 
nary SMBH systems. 

We would like to thank Philipp Kronberg for providing the 
radio map of 3C9. Support for this work was provided by the 
NSF through grant AST-007 1181. 

A. MORE ON NON-CIRCULARLY SYMMETRIC POTENTIALS 

Many observed gravitational lenses do not produce simply 
two bright images as a circularly symmetric potential would 
predict. Instead, four or more bright images are often ob- 
served. Two effective and common explanations are: 1) that 
the lensing galaxies are indeed not circularly symmetric but 
rather elliptical, and 2) that the lensing galaxies lie in an ex- 
ternal gravitational shear. In this section we provide a cursory 
analysis of how ellipticity and shear will modify the results 
of the previous sections. In addition, we treat the case of an 
off-center SMBH in the large source impact parameter limit. 

We have seen in the case of circularly symmetric potentials 
that, for large impact parameters ((3), the magnification of im- 
age D is essentially independent of the lensing properties of 
the galaxy that contains the image. We will show that neither 
ellipticity nor shear, will significantly alter the magnification 
of this image or the primary image, and therefore, the opti- 
cal depths and statistics pertaining to observing these images 
calculated for the circularly symmetric potentials do not need 
to be altered. The properties of any other faint images, how- 
ever, will be complicated by the addition of the more complex 
potentials and our previous calculations may not necessarily 
hold. Similarly, an off-center SMBH is shown not to alter the 
magnification of image D in the large source impact param- 
eter limit, but other configurations and other images may be 
affected. 

A.l. Ellipticity 

It has been demonstrated (Evan s & Hunterir2002l) that an 
elliptical lensing potential derived from a cored isothermal 
sphere increases the maximum number of images produced 
from three to five. For our PM+CIS potential, considering 
ellipticity does not actually increase the maximum number of 
images, but r ather com plicates the caustics in the source plane 
considerably (Mao et al. 2001). We expect that the magnifica- 
tion of the SMBH image, where it exists, and of the primary 
image will be affected little by this change, but to demonstrate 
this quantitatively we must return to the lens equation. Recall- 
ing from Equation^ but this time in two dimensions, the lens 
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equation is: 



/3 = 0-V*. 



(A-l) 



In order to describe a galaxy with an elliptical density distri- 
bution, we will modify the circularly symmetric cored isother- 
mal sphere lensing potential to (Naravan & Bartelmann 1996, 
Section 3.5.1): 



*elli P tical(O x ,ey) = S is ^£ + (1 - e)0 2 + (1 + e)02 (A-2) 

In this equation the variables t and 9 y are the angular dis- 
tances along the principal axes of the ellipse, and e is the el- 
lipticity of the potential. The cored elliptical lensing potential 
is not entirely realistic since, for large ellipticity, it results in 
dumb-bell shaped mass-density contours. For small elliptic- 
ity, however, the contours are generally elliptical and the po- 
tential is quite acceptable. Combining this with the point mass 
potential yields the complete lensing potential for an elliptical 
galaxy with a SMBH: 



n0,A) = 0^Jel + (\-e)9l + (\ + e)9l+9l M \nJel + e 1 y . 



(A-3) 



Using the lens equation, we find for the components of f3: 



0x = X 



Py = 0y 




9 2 c + (l-e)9 2 x + (l + e)9 2 * + y 



(A-4) 



(A-5) 



Since we are interested only in magnification of images for 
the present purposes, we can forgo completely solving the 
lens equation. For non-circularly symmetric potentials, the 
magnification of any image is given by ( Scheider et al. 1992, 
Equation 5.21): 



fi =det.4, 



(A-6) 



where A is the Jacobian matrix which describes the lens prop- 
erties, and is: 



• A " de 



d x (3 x d y (3 x 

O x Py dy(3 y 



(A-7) 



Unfortunately, writing out A and calculating its determinant 
is not particularly illuminating as it is difficult to see past the 
algebra. As a simplification, we consider the limiting cases 
of 9 <C (6pm < 6>c) and the opposite extreme 9 ^> (9 PM < 0c). 
The former happens to correspond to the region of the im- 
age plane where SMBH images are expected to be observed 
because, from our earlier analysis, we might reasonably as- 
sume c /0sis - 10 _1 , PM /0sis ^ 10" 2 , and for SMBH im- 
ages, 0/0sis — 10~ 3 . The latter corresponds to the region 
where the primary image is observed since 0/0sis ^ 1 for such 
images and the other parameters, of course, do not change. 

For the case of the inner region of the image plane, where 
the SMBH image is expected to form, the components of (3 
reduce to: 



lim (3 X = X 

#-C#PM 



lim j3y = Oy 

0*C#PM 



r al -i 

"PM 


[ oi+e 2 y\ 

r (ft i 

P PM 


02 + 02 _ 



(A-8) 
(A-9) 



and thus the magnification for an image in this region is given 
by: 



I^smbh ~ det.4 '■ 



9 4 
7 PM 



9 4 
7 PM 



02 + 02)' 



■f- (A-10) 



The magnification of these images would not, therefore, be af- 
fected by the ellipticity of the lens galaxy. It should be noted, 
however, that this approximation breaks down in this region 
if 0c/0sis is sufficiently small (~ 2 /0p M ). 

For the case of the outer region of the image plane, where a 
primary image forms, the components of (3 reduce to another 
pair of relatively simple expressions: 



Q-»e c 



lim (3 y = 9y 

0»0c ' 



7SIS 



(1-6) 



(l-e)0|+(l + e)02 

0sisd+e) 
(l- e )02 + (l+e)02_ 



(A-ll) 



(A-12) 



Along the axes, 9 X and 6 y , the magnification in this region is 
given by: 



M Primary ~ UCtA ; 



1- W when0,. = O, 



sd-e) 



when 9 X = 0. 



(A-13) 



We can consider that, in effect, the addition of ellipticity sim- 
ply modifies the value of 0sis depending on the position of the 
primary image. Since the value of e is at the very greatest 1, 
and is usually much smaller, we can ignore this contribution 
as a minor correction. 



A.2. Gravitational Shear 

The lensing galaxy might be residing, on the other hand, 
in the presenc e of an external g ravitational shear. It has been 
demonstrated (Evans & Hunter 2002) that the addition of an 
external shear to a cored isothermal sphere potential increases 
the maximum number of images produced from three to five, 
just as did the consideration of ellipticity. For our PM+CIS 
potential, the maximum number of images increases again 
to six, however four images is a more probable outcome for 
most configurations, as shown in Figure[21 Again we expect 
that the SMBH image and the primary image will be affected 
little by this change. Returning to the two-dimensional lens 
equati on, we begin by ch aracterizing an external shear with 
JNaravan & Bartelmannlf9 96. Section 3.5.2): 



*shear(0„0v)=^(0?-0v)- 



(A- 14) 



Here 7 is the magnitude of the shear and, as in the case of 
the ellipticity, the angular distances 9 X and 9 y measure along 
the axes of the shear. In order to consider how this affects the 
images, we begin by adding this contribution to the potential 
for a PM+CIS giving: 



*(0„0,) = 0SIS 



0? + 02 + 02 + 6^ In 



W 



V^) + \i8l- 



9;). 



(A-15) 
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From the lens equation we again can calculate the components 

of/3: 



&=0v 1 



(3y = 6 y \ 1- 



0SIS 



9 2 

7 PM 



ei+e 2 y +ei % +e y 



7 SIS 



9 2 

7 PM 



0? + 2 + 2 



61+ & 



-7 (A-16) 



+ 7 (A-17) 



and for convenience we will define an angle ui which mea- 
sures the position of the image with respect to the positive 9 X 
axis such that: 

$ x = 0cosui (A- 18) 

y = 0sinw. (A- 19) 

Calculating the determinant of the Jacobian matrix A and 
grouping the terms by powers of 7 yields, after some sim- 
plification, 



det„4 = 



1- 



9 4 

7 PM 



q2 a2 
'SIS^C 



, + fl z 



26l + t 



+ 7Cos2cj 



/)2 

2 ' 



7 SIS" 



7SISC 



:+ 2 ) 



3/2 



,3/2 



"7 



(A-20) 



Following the same approximations here as we did for the 
elliptical case, we consider the two regions of the image plane 
where < (0 PM < 0c) and > (0 PM < 0c)- For the inner 
region of the image plane where an SMBH image is expected 
to be observed, and -C (0 P m < 0c), we substitute, as we did 
before, the values of 0c/0sis — 10 _I , P m/< 7 'sis — 10~ 2 , and 
10~ 3 . In this limit, the above expression reduces to: 



7 SIS 



lim (det„4) 



9 4 

7 PM 



a 

PM 



& 4 +2 7 cos2^ g2 



"7 



(A-21) 



Typical values of the shear parameter are generally of order 
7 < 0.1 (|Casertano et al. 2003). Thus, it is sufficient to keep 
only the leading term of the expression for the magnification: 



^SMBH ' 



M 
^PM 

4 ■ 



(A-22) 



This is identical to the elliptical result for this region and we 
conclude that shear will not substantially alter the magnifica- 
tion of an SMBH image for such a lens configuration. 

For the case of the outer region of the image plane, where 
a primary image forms, and ^> (0 P m < 0c), we proceed by 
repeating the above simplification, but with typical values of 
1 . The determinate of A is then reduced to: 



7 SIS 



7SIS 



lim (det„4) — ► 1 — + 7cos2o; 

0»0 C 

and using 7 < 0.1, 



7 sis 



"7 



^Pr, 



7 SIS 



1 — [l-7Cos2w] . 



(A-23) 



(A-24) 



The angular dependence in this equation makes sense, be- 
cause symmetry dictates that two points separated by 180 de- 
grees must have the same magnification. As was the case for 
the elliptical potential, this can be considered to be a small 
correction to 0sis and, therefore, the presence of an exter- 
nal shear has not altered the ratio of the magnification of a 



SMBH image to a primary image in these limiting cases. Fig- 
ure^]plots the image magnifications for a representative sys- 
tem with 7 = 0.1 as function of source position for a path ex- 
tending away from the center of the lens. Comparison with 
the top-left plot in Figure [5] which has no shear but is other- 
wise the same configuration, demonstrates the minimal effect 
of shear on image magnifications for both core images in this 
case. 

A3. Off-Center SMBH 

Anothe r possi ble departure from circular symmetry, raised 
in Section R-,3.31 is that the SMBH may not be located in the 
center of the host galaxy. Returning to the lensing potential 
given in Equation|6j we can explore the consequences of such 
an occurrence by inserting an offset, 9q, between the point 
mass and the galaxy potential, 



y = 9 sl5 ^\e\2+e 2 c +el M ]n\e-0 o 

This yields the lens equation, 
9 



[3 = 9-9: 



sis 



9 2 \ + 9 2 c 



,2 Q-Oq 

PM |6/-6>o| 2 ' 



(A-25) 



(A-26) 



Following a similar approximation here as we did for the el- 
liptical and shear cases, we consider the region of the image 
plane which corresponds to large source impact parameters. 
In this limit, we expect that an image will form near the point 
mass, and thus \8 — 0q\ will be small. The last term in the 
lens equation, therefore, will dominate and we can simplify 
the lens equation to, 



/3«-0 



e-Qo 

9-9 P- 



and rearrange to solve for the image position, 

/3 



9 



SMBH ■ 



00-^Mpi- 



(A-27) 



(A-28) 



The position of the SMBH image as a function of source im- 
pact parameter is merely offset from the original location by 
the same amount as the SMBH itself. The magnification of 
the image, 



(J'SMBH 



9 4 

7 PM 



(A-29) 



will also have the familiar form, but with respect to the new 
SMBH location. 

A. 4. Remarks 

Neither ellipticity nor external shear (nor, under proper cir- 
cumstances, an off-center SMBH) significantly change the 
probability of observing lensing by a SMBH. The cross- 
sections, optical depths, and statistics we calculated for the 
circularly symmetric models, therefore, remain a reasonably 
good representation even under these more general circum- 
stances. Other core images, however, do not necessarily lie 
in the limiting case regions of the image plane and may be 
affected by ellipticity and shear. We cannot assume that the 
results of the circularly symmetric calculations can be applied 
to deduce their statistics for these configurations, but postpone 
further analysis for future work. 
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FIG. 1 . — Positions and magnifications of the images produced by PM+CIS model (top) compared to the cored isothermal sphere (bottom). Image positions and 
the source position are given in units of the singular isothermal sphere Einstein ring radius, $sis- For this figure we chose #c/$sis = 0.01 and #pm/#sis = 0.03. 
The primary image is labelled "A" and the secondary images are labelled "B," "C," and in the case with the point mass, also "D." The dashed lines indicate the 
properties of the secondary image of an isolated point mass. 
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FIG. 2. — Einstein radii ratio #pm/#sis is plotted as a function of source redshift for a lens galaxy at zl = 0.5 (top), and as a function of lens redshift for a source 
galaxy at zs = 2.0 (bottom). The empirical relationship between SMBH mass and host galaxy velocity dispersion is used and three values of velocity dispersion 
are shown. A standard cosmology of Qm = 0.3, C2a = 0.7, and Hq = 70km/s/Mpc is used. 




FIG. 3. — Physical distance, r c , corresponding to the angular core radius ratio, #c/6sis = 0-1, is plotted as a function of source redshift for a lens galaxy at 
Zl = 0.5 (top), and as a function of lens redshift for a source at zs = 2.0 (bottom). Other values of the ratio can easily be deduced because the curves scale linearly 
with #c/6sis- Three values of velocity dispersion are shown. A standard cosmology of Qm = 0.3, Qa = 0.7, and Ho = 70km/s/Mpc is used. 
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FIG. 4. — Series of plots indicating the influence of PM+POW parameters on lensing. The left column is a series of plots illustrating the progression, from 
top to bottom, of the magnification curves for the PM+POW potential as the sharpness of the break, a, is increased while holding constant r\ = 1, 7 = 0.1, 
^b/^sis =0.01, and #pm/#sis =0.03. The profile of the normalized surface brightness, Y,powlh> is drawn along side for reference. On the right, the progression 
illustrates increasing the inner power law index, 7, while holding constant a = 2, r\ = 1, #b/$sis = 0-01, and Opyi/dsis = 0.03. The final frame, with 7 = 1, is 
effectively a singular isothermal sphere plus SMBH. 
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FIG. 5. — Series of plots illustrating the progression, from top to bottom, of the magnification curves as the core radius is reduced. The left column depicts 
the progression for the PM+CIS potential. The right column depicts the progression for the PM+POW potential where a = 2, rj = 1, and 7 = 0.1. For both cases, 
^pm/^sis = 0.03. The behavior is qualitatively similar in both cases. As the core radius is reduced, the region of four images moves outward from the lens galaxy 
and reduces in area until it disappears at a finite core radius (not shown). 
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FIG. 6. — Normalized cross-sectional area, A/(7r6> 2 IS ), of the anulus in the source plane where four images are produced is plotted as a function of core radius. 
The left graph is for PM+CIS, the right is for PM+POW. The curves represent different values of 9pm/^sis- Beginning with the lowest curve in each graph and 
working up, #pm/0sis = { 10 _1 , 10~ 1,2 , 10~ 14 ...} . The solid line marking the upper bound of all curves represents the area of the corresponding three-image 
region for a cored isothermal sphere with no SMBH. The reduction of the curves from the solid line, particularly for small core radii, represents core image 
"swallowing" due to the SMBH. 
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FIG. 7. — Normalized cross-sectional area, A/(rr6| IS ), in the source plane is plotted as a function of core radius for 8pu/9sis = 0.03. The curves represent 
different values of the flux ratio, R, ranging from 10 8 down to 10°. The left column of graphs is for PM+CIS, the right is for PM+POW. The first row is the 
cross-section for two images, the second row is the cross-section for three images, and the third row is the cross-section for four images. 
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FIG. 8. — Normalized cross-sectional area, A/(7r#| IS ), in the source plane is plotted as a function of core radius for a dynamic range, R = 10 6 . The curves 
represent different values of $pm/#sis ranging from 10' down to 10 -4 . The left column of graphs is for PM+CIS, the right is for PM+POW. The first row is the 
cross-section for two images, the second row is the cross-section for 3 images, and the third row is the cross-section for four images. 
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FIG. 9. — Optical depth, r, for lensing of a source at redshift, zs = 1. 1, expressed as a function of core radius. The curves represent different values of the flux 
ratio, R, between the brightest image and image of interest. The left column of graphs is for PM+CIS, the right is for PM+POW. The first row is T5, the optical 
depth for observing a lens with at least two images. The second row is t-j, the optical depth for at least three images. The third row is T4, the optical depth for 
four images. Note that for T3, the curves above R = 10 5 saturate and overlap, with R = 10 3 5 being the first to significantly break apart. The same is true for T4, 
however in this case, the curves above R = If) 6 saturate and R = 10 5,5 is the first to part significantly. The saturation is the result of the finite area of the region 
between the caustic circles in the source plane that produces more than two images. If the source redshift is increased to zs = 2.0, then the optical depth curve for 
a given R is increased by approximately a factor of 3. A standard cosmology of Qm = 0.3, Q\ = 0.7, and Hq = 70km/s/Mpc is used. 
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FIG. 10. — Number of lenses predicted to be observable in the FIRST radio source catalog. The top row is for a search with limiting sensitivity S = 1/iJy 
and is plotted as a function of core radius. The curves are given for the number of lenses with at least two, three, or four images detectable (N2, N3, and N4, 
respectively). A 'outside ls tne number of lenses with two images detectable where the second image would not be present without the SMBH. In the middle row, the 
plots represent the fraction of observable lenses with detectable third images. The curves, in this case, represent different limiting sensitivities from S = 1000/Ujy 
down to S = O.OljuJy. The bottom row of plots is, similarly, the fraction of observable lenses with detectable 4th images. For all rows, the left column of graphs 
is for PM+CIS, the right is for PM+POW. The vertical dashed lines are drawn for the values of 8q used for the discussion in the text. 
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FIG. 1 1. — Contour plot of the radio source 3C9. The total flux is approximately 190 mJy. The location of two nearby galaxies are also indicated. Galaxy 1 is 
about 10 arcsec away from the source, and galaxy 2 is about 3 arcsec away. 
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FIG. 12. — Simulated lens image of 3C9 produced by a 10'Mq SMBH at the center of galaxy 1. The position of the black hole in this plot is the upper left 
corner. The total flux of the image is approximately 0.03 njy. 
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FIG. 13. — Simulated lens image of 3C9 produced by a 10 9 Mq SMBH at the center of galaxy 2. The position of the black hole is marked by the cross. The 
total flux of the image is approximately 3 nJy. 
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FIG. 14. — Lensed images of 3C9 produced by a binary SMBH system in galaxy 2. Panel A shows the image produced by a single 10'Mq SMBH at the 
position of the cross. Panels B, C, and D are the images produced by two 5 X 10 Mq SMBHs separated by 0.010pm (2.5 pc), 0.020pm (5.0 pc), and 0.03#pm (7.5 
pc), respectively. The total flux for each panel, from A to D, is: 3.1 nJy, 4.0 nJy, 3.5 nJy, and 2.3 nJy. 
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FIG . 15. — Stylized schematic of the caustics in the source plane for the sheared PM+CIS potential. The number of images produced in each region is given. 
The relative sizes of the caustic curves are dependent on the choice of parameters and the region producing six images is typically not present for small core radii. 
The dashed line indicates the approximate path of source positions for which the image magnifications are plotted in the bottom panel. In the bottom plot, the 
source position is given in units of the singular isothermal sphere Einstein radius, $sis, and we use #c/#sis = 0.1, Spm/^sis = 0.03, and 7 = 0.1. 



